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Abstract 

We characterize the existence of the solutions of the truncated 
moments problem in several real variables on unbounded supports by 
the existence of the maximum of certain concave Lagrangian functions. 
A natural regularity assumption on the support is required. 
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1 Introduction 

The present paper is concerned with the truncated problem of moments in 
several real variables, in the following context. Let n G N and fix a closed 
subset T 7^ of M", a finite subset / C with G / and a set = {gi)i^i 

of real numbers with go = 1, where Z+ = N U {0}. Typically a problem 
of moments [1] requires to establish if there exist Borel measures z/ > on 
M", supported on T, such that |(iz/(t) < oo and jr^fdvif) = gi for all 
i & I. As usual f = fi ■ ■ - f^n where t = {ti, . . . ,tn) is the variable in R" and 
i = {ii, . . . ,in) is a multiindex. In this case we call u a representing measure 
of g, and gi the moments of u. We are interested in those measures u = fdt 
that are absolutely continuous with respect to the n-dimensional Lebesgue 
measure dt = dti ■ ■ ■ dtn, in which case we call / a representing density of g. 
Namely the (class of equivalence of the) Lebesgue integrable function / is 
> almost everywhere (a.e.) on T, has finite moments of orders i E I and 



[ ffit)dt = g, {lel). (1) 
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Given partial information in the integral form f f pdt = g.i about rep- 
resenting densities f on a probability space (T,pdt), endowed with a refer- 
ence density p, does not determine them uniquely. An approach favorite to 
physicists and statisticians is to choose that particular density f*, minimiz- 
ing the entropy functional h{f) = Jj,{fliai)pdt amongst all solutions of the 
moments constraints. This uniquely selects the unbiased probability distri- 
bution f* (that proves to have the form f^,(t) = e^»6-f^»*') on the knowledge 
of the prescribed average values gi of f, where t is considered as a T- valued 
random variable with repartition p [6], [9], [18], [20]. Under suitable hy- 
potheses, f* turns to exist, even for measures more general than pdt. A 
main tool to this aim is Fenchel duality [S], [21], [2S], [27], that deals with 
minimizing convex functions h : X ^ MU{cxd} on convex subsets of lo- 
cally convex spaces X, in connection with the dual problem of maximizing 
— where h* : X* — > MU{cxd} is the convex conjugate of /;., called also 
its Legendre- Fenchel transform [26], [27], defined on the dual X* of X by 
h*{y) = s\iY){{x,y)—h{x):h{x)<oo}. Typically inf /i = max(— /i*) and, briefly 
speaking, minimizing J^f Infpcit as above is to find A* = (A*)jg/ maximizing 
— Ylii&iQi^i ~ Jj.e^iei^i^'' pdt. Many results exist in this direction [3], 
[3] - [H], [IS], [H], [21] - [23]- Additional hypotheses are always necessary 
when the conclusion inf h = min h is sought for, since there are g for which 
the primal attainment fails [16], [17] although problem ([1]) has solutions. 

By Theorem |3] we prove that the feasibility of problem ([1]) is equivalent 
to the boundedness from above sup L < oo with attainment sup L = max L 
for the concave Lagrangian function L. This holds no matter whether inf h 
is attained or not (the general theory still provides us with inf h = maxL). 

Initiated by Stieltjes, Hausdorff, Hamburger and Riesz, the area of the 
truncated problems of moments nowadays knows various other approaches, 
based for instance on operator methods or sums-of-squares representations for 
positive polynomials [TU] - [H], [IH], [25]- Although important, these topics 
remain beyond the aim of this work, focused on our mentioned Theorem [3J 

The author got the idea to consider L instead of h from the works [5] 
where a similar characterization exists, and [16], [17], drawn to his attention 
by professor Mihai Putinar. Our statement and proof are rather general, 
independent of these cited works. 
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2 Main results 



We remind that a linear Riesz functional ip^ [12] associated to a set 7 = 
of real numbers 7^ for J C is defined on the polynomials p from 
the linear span of . . . X^" where i = (ii, . . . , i„) G J by v^-yX* = 7j. One 
calls if^ T-positive [12] if ip^p > whenever p{t) > for all t ^ T. If 7 has 
representing measures z/ > on T, 9?^ is T-positive since ip^p = Jj.pdu for 
any such polynomial p. In the full case J = Z" the T-positivity condition 
is sufficient for the existence of the representing measures, by the Riesz- 
Haviland theorem [15]. An analogue of this theorem [12] for the truncated 
case I = {i: \i\ <2k} characterizes the existence of the representing measures 
by the existence of T-positive extensions of (p^ to the space of polynomials 
of degree < 2k + 2. For later use, we state below a version of these results 
(Theorem [1]) and a Fenchel theoretic result of dual attainment (Theorem |2]). 



Definitions We call T regular P] if for any t & T and e > the Lebesgue 
measure of the set {x G T : ||x — t\\ < e} is positive. As usual = 
(Z]r=i ^fy^'^- ^'^^ ^ G set (Tj = {j G Z" : jfc = either or z^, 1 < /c < n}. 
We call / regular [1] if C / for all i G /. Define T, G CR^ {N = card/) 
by r = {7 = (7j)ie7 : 3 measures z/ > on T with Jr^fduit) =7^, i E 1} and 
G = {7 = (7,)ie/ ^ : 3 / G L\{T,dt) such that j^tf {t)dt = 7i, i G /}. 
The notation LP{T,fi), L^in) for measure on T, 1 <p< 00 has the usual 
meaning. In particular L\{T, fj,) is the set of all / G L^{T, fi), / > /i-a.e. For 
7 = (7j)je/; f-y is ^f^^ linear functional defined on the span Pj C . . . , X„] 

of all X' with zG/ by (^^X* = 7i. Set e, = (0, . . . , 1, . . . , 0) for l<L<n. 

By [Theorem 6,[1]] the convex cone G is the dense interior of the cone F. 

Theorem 1 [Theorem 7,[1]] Let T C 6e a closed regular set, I C Z" a 
/imie regular set and g = {gi)i^i a set of numbers with go = 1. Then g & G 
(PgP > for every j9 G P/ \ {0} such that p(t) > for all t E T. 

Theorem 2 [Corollary 2.6, [8j] Let T he a space with finite measure fi>0, 
1 < p < 00 and ai G L'^^jj), (?i G M for i E I = finite where ^ -f- ^ = 1. Let 
: M— 7- (—00, 00] be proper, convex, lower semicontinuous with 0|(o,oo) < 00. 
If there arexEL'^^jj), x>0 a.e. such that (pox E L^i^jj) and J^aiX dfi = gi, 
then the quantities 

P = mi{J (f){x{t))d^{t) : x E L^(/i), x > a.e., (j)ox eL^(h),J aixd^ = gi\/i }, 
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are equal, —oo<P = D<oo and the maximum D is attained. 

Theorem [3] is a reminiscent to [Theorem 4, [3]], where J^flnfprft is 
minimized subject to jj,ffpdt = Qi under stronger hypotheses on p, like 
pit) ~ e~^ll*ll'' with p > 2k (to fit the notation in [3], let a = 1 and our 
f := pi, whence Lp^a,g{><) = L{X- Xo) + l, with Aq = (Aoi)ie/ where Aoi = 5i,o 
and 6ij is Kronecker's symbol, 6ij = 1 if i = j and if i^ j). Although 
we do not obtain here the existence of a maximum entropy solution f^, our 
present hypothesys on p are weaker, while condition g & G still characterized 
in Lagrangian terms. Our proof below relies on Theorem [1] ( [Theorem 7,[1]]) 
and Theorem [2] ([Corollary 2.6, [8j]). 

Theorem 3 Let T C M" be a closed regular set. Let I C be a finite 
regular set such that maxi^j \i\=2k where k&N. Assume 2ke^&I {l<L<n). 
Let g = {gi)i(^i be a set of numbers with g^ = 1. Fix p G L^{T, dt), p > a.e. 
The following statements (a) and (b) are equivalent: 

(a) There exist functions f G L^{T,dt) such that jr^ \f\f{t)dt < oo and 

jj'f{t)dt = g, {zeiy, 

(b) The functional L : — 5- M U {— oo} defined by 




'p{t)dt, A = (Ai)ie/ 



is bounded from above and supL is attained in a (unique) point A*. 

Proof Since L(0) > — oo, L ^ — oo. Since go = 1, each of the conditions 
(a) and (b) implies that T has positive Lebesgue measure, finite or not. Hence 
by means of Jensen's inequality one can show that L is strictly concave. Then 
whenever supL is finite and attained at some point A*, this A* is unique. 

(a) =^ (b) The regularity condition on T is not necessary for this implica- 
tion. Let /i = pdt be the measure on T with density p := pe~^'=i*' . Then 
< p(T) < oo. Since ([1]) has a solution /, then f := f / p satisfies 

/ ff{t)dpit) = g, (tel). (2) 
Jt 
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By [Theorem 2.9, [H]], see also [Lemma 4, [1]] for (3 = 0, problem ([2]) has 
also a solution /o G L°°(T) with /o > a.e. The conclusion supL < oo 
may hold either directly by Theorem [21 or by an elementary argument as 
shown below. Let x = foit) a.e. and y = ||/o||oo + 1 in the inequalities 
—e~^ < xlnx < ylny for < x < y, y > 1, then integrate with respect to fi. 
Hence /oln/o G L^{T,n). Fix A = (Ai)jg/. Let x = fo{t) and y = Y.i&i^i'^' 
in the simple version xhix — x > xy — oi Fenchel's inequality , then 
integrate. It follows, using ([2]) for /o, that 

/ /oln/orf/i- / /oc?/i > V^iA, - / e^'e^^'*'c//i(t) = L(A-Ao) + y'^iAoi 

where Aq = (Aoi)iG/ with Aqj = TTi=i^i,'ikc, and 5ij is Kronecker's symbol. 
Since A was arbitrary, we get sup;)^L(A) < oo. Now for the attainment 
supL = maxL, we need Theorem [2] as follows. Use \tj\ < (YT^^itf^Y^'^'^ , 

n n 

|f I = . . . < (^tf ^ ^yi\/2k < ^^2fe ^ ^ ^1 -J < 2k) 

and i/+l<e^ for y = YJl=i'tf to get |f |c?/i(t) < J^pdt <oo for ieJ. Then 
let: T = T, the measure /i = ptit, p = oo, the moment functions ai{t) = f 
and the integrand be defined by = xlnx for x > 0, 0(0) = and 
0(x) = +00 for X < 0. The feasibility hypotheses is fulfilled by x = /q. The 
convex conjugate 0*(|/) = snp^yQ^xy — xlnx) of is given by 0*(?/) = e^~^ 
for y e M. We get the attainment D = sup £ for £(A) =L(A— AQ)+^.gj SfiAgj 
where Aq = (Agjig/ with Aq^ = Aoi + 5i,o- Thus we obtain a A* such that 
supL = L(A*). 

(b) ^ (a) Let A* G such that supL = L{X*). We prove that (pg 
satisfies the positivity condition in Theorem [TJ Let p = X^ie/ P ^ 

be arbitrary such that p(t) < for t G T. The vector A := (Aj)ig/ is then 
^ 0. For any r > 0, set e,.(t) = e^^^ei^'*\ Thus e,.(t) < 1 for t G T. 
Then the integral term fj,erpdt of L(rA) = rj^ieidi'^i ~ frpCrpdi remains 
bounded as r — )• oo. Hence ipgp = ^jg/fi'jAj < 0, for otherwise the linear 
term npgp of L(rA) would give sup L = oo that is false. Assume that (fgp = 0. 
Then the restriction of the function L to the half-line i := {rA : r > 0} is 
given by the function r — J^, Crpdt. This function is finite, bounded and 
strictly monotonically increasing on (0, oo). Use to this aim that < < 1, 
Jrp pdt < oo, Cr = e^^ with p < and L\i is strictly concave. Then a finite 
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limit lim L{rX) = sup L exists, in particular sup |L(rA)| < oo. For a>0, 

r— >-oo ^ r>l 

oo>L(A* + aA) = ^^iA* + a^^iAi- f e^'^' ^''^ e''^^<^' ^^'^ p{t)dt 
iei i&i 

> y'^iA* + ^ - - / e^^e' p{t)dt = L{X*) = maxL > L(0) > -oo 

because 9i^i ~ ^ Sie/ -^i^* — ^ t E T. Hence L is finite on 

every point of the half-line {A* + aX}a>o- Note that A* cannot be colinear 
with A due the behaviour of L on i: firstly, A* ^ £ because L reaches its 
global maximum only in A* while L\£ increases strictly along £ as r — )■ oo. 
Also A* ^ {0} U (— £), for otherwise the concavity of the restriction L\^x '■ 
MA — i- {— oojUMofLto the line MA would imply, for some r > with 
A* = -rA, that L(rA) > L(0) = i:(|(A* + rA)) > i(L(A*) + L{rX)), whence 
L(A*) < L(rA) < supL|^ < supL = L{X*) that is impossible. Thus A* ^ MA. 
Then a 2-dimensional drawing shows that for every r > 1 there is a unique 
point Xr of intersection of the segments (A*,rA) and (A, A* + A). Write to 
this aim Xr = sX* + (1 — s)rX = s'X + (1 — s')(A* + A) with coefficients 
s = Sr, s' = s'j., use the linear independence of A*, A and get s = (r — l)/r, 
s' = 1 — s whence s, s' G (0, 1) and lim^_j.oo s'^ = 0. Then limr_^oo Xr = X* + X. 
The concavity (and hence, continuity [27]) of L on the segment (A, A* + A] 
gives limr^oo L{xr) = L{X* + A) < ^^(A*) with strict inequality, because the 
point A* of maximum of L is unique. But L{xr) = L{sX* + (1 — s)rX) > 
sL{X*) + (1 — s)L(rA) and letting r — )■ oo we derive, using limr_).oo Sr = 1 and 
sup^>]^ |L(rA)| < oo, that limr^.oo -^(a^r) > L{X*). We got a contradiction. 
Then (fgp < 0. The feasibility of problem ([1]) follows then by Theorem [H □ 

Remarks Since A* may be on the boundary of domL := {A : L{X) > -oo}, 
one cannot prove (b) =^ (a) by derivating under the integral in A*, and the 
/i- minimization may fail [17]. Additional hypotheses may compel A* to be 
interior to dom L [16] in which case the entropy minimization can be obtained 
[21], providing the particular solution =e^'6^'^**\ see for instance [3]. 
For example let T = M", I = {i : \i\ < 2k} and p{t) = e~"*" . By Theorem 
ini problem ([1]) is feasible if and only if L is bounded from above and attains 
its maximum in a point A*, even when a minimum entropy solution does not 
exist. By Fatou's lemma and Lebesgue's dominated convergence theorem, 
/o := e^i»i^2fc '"'i* ]2as finite moments of order < 2k, we can get J ff()dt = gi 



6 



for \i\ < 2k and / tf^ fodt < g2kc, (1 < < ''^), but the equalities ([T]) may 
fail for 1^1 = 2k ^7\. By integration in polar coordinates, the homogeneous 
polynomial p := '^\i\=2k ^i-^^ shown to always satisfy p(t) < on R"; if 
moreover p{t) < for all t^O, then A* is interior to domL and /o is indeed 
a solution of problem (II]), /o = /*. We omit the details and refer the reader 

to m, m- 

Note also that whenever p is at our disposal, various choices may be tried 
[3] to facilitate the numerical maximization of L = Lp. 
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